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ABSTRACT 


This  study  treats  the  ducted  propeller  with  finite 
blade  nuniber  at  zero  angle  of  attack  in  a  uniform,  incom¬ 
pressible,  inviscid  flow.  The  approximations  of  a  lightly 
loaded  propeller  and  of  thin  airfoil  theory  are  made.  In 
the  absence  of  thickness  effects,  appropriate  vortex  dis¬ 
tributions  represent  the  blades,  the  shroud  and  their  res¬ 
pective  shed  vortices.  By  means  of  Fourier  2Uialysis  of  the 
velocity  field  in  propeller  fixed  coordinates,  the  problem 
for  an  arbitrary,  radial  blade  circulation  distribution  is 
reduced  to  a  form  similar  to  the  ring  wing  integral  equations 
of  j.  weissinger.  The  kernels  are  not  the  same  for  the  two 
cases  except  for  the  zeroth  harmonic.  The  equation  for 
this  harmonic,  which  corresponds  to  a  generalized  actuator 
disk  formulation,  is  identical  to  that  of  an  equivalent, 
axisymmetric  ring  wing.  The  effect  of  blade  number,  blade 
circulation  profile  and  strength,  propeller  advance  ratio 
and  tip  clearance,  and  location  of  the  propeller  plane  are 
found  and  may  be  evaluated  over  a  range  of  parameters  from 
tables  provided.  A  numerical  example  is  given. 


ii 


TABLE  OF  CONTENTS 


INTRODUCTION  1 

CHAPTER  (BIE  -  Basic  Fomulatlon 

1.1  Physical  Problem  and  General 

Background  3 

1.2  Coordinate  System  3 

1.3  The  Mathematical  Model  8 

1.4  Velocity  Influence  Functions  13 

1.5  Calculation  of  the  Induced  Velocities  16 

1.6  Formulation  of  Basic  Equation  for  y  19 

CHAPTER  TWO  -  Derivation  of  Shroud  Vortex  Distribution 

2.1  Periodic  Expansion  of  Shroud  vortex 

Distribution  i.? 

2.2  Reduction  of  Shroud  Contributions  2: 

2.3  Decomposition  of  Governing  Equation 

for  7  28 

2.H  Splitting  of  Governing  Conqplex  Equation  32 

2.3  Discussion  of  Coupled  Equations  36 

2.6  Shroud  and  Propeller  Loading  40 

CHAPTER  THREE  -  Solutions  for  Duct  Leading 

3.1  Fundamental  Solution  43 

3.2  Effect  of  Parameters  or  the  Fundamental 
Load 

3.3  Limiting  Case  of  Infinite  Advance  Ratio  6l 

3.4  Infinite  Blade  Number  64 

ill 


CBAFTEK  THREE 


3 >5  G«n«ral  solution  for  Higher  Harmonica 
3.6  Exaa^le 


REFERENCES 


•  •  • 


0„(«2) 


ai4>lltudes  of  sine  tezn  of  Fourier  harnonlc 

amplitudes  of  sosine  term  of  Fourier  harmonic 

amplitudes  of  coo^lex  Fourier  hansonic 

distance  between  vortex  elesMnt  and  field 
point 

f;inctions  of  complex  kernel  of  intermediate 
form  of  governing  equations 

Riegels  function 

vector  influence  functions  for  axial,  radial 
and  circumferential  vortex  eleewnts  respectively 


J 

V 


lntegi-«ia  of  choxdwlse  vortex  distribution 
and  decoupled  kernel 

imaginary  part  of  complex  function 

propeller  advance  ratio, 

kernel  of  final  governing  equations 

blade  nusOser 


P 


"forcing  function"  of  higher  propeller 
harmciiice 


Legendre  function  of  second  kind  and  half  order 

B  radius  of  the  shroud  camber  surface  in  the 

propeller  plane  used  as  reference  length 

real  part  of  complex  function 


R  propeller  semi -diameter 

P 

®n  ^°n+l/2'^®n-3/2^ 

’’^n  ^®n+l/2"°n-3/2^ 

U  uniform  free  stream 


V 


c 

shroud  chord 

^saN'^toN 

fvmctlons  dependent  on  higher  harmonics  of 
shed  vorticity 

i.i.K 

unit  vectors  in  the  cartesian  syston  (x.y.z) 

unit  vectors  in  the  cylindrical  system  (x.r.Q) 

i 

^1 

j 

dummy  summation  index 

i 

blade  index  number 

a 

rank  of  Fourier  harmonic 

n 

order  of  Legendre  function 

P 

static  pressure 

2 

fluid  velocity 

t 

time  variable 

X.X.& 

cylindrical  coordinates  fi.  ed  in  propeller 

S  3  S 

cylindrical  coordinates  of  ,  .int  on  shroud 
surface 

x.y.z 

Cartesian  coordinates  fixed  in  th>'  propeller 

s  s  s 

Cartesian  coordinates  of  point  on  shreud 
surface 

X.  ... 

nondimens ional  coordinates  v;ith  respecc  to 

R  ;  e.g.,  X  =  x/R 

X,  ... 

nondimens ional  coordinates  with  respect  to 

Rp  ;  e.g.  X  --  x/Rp 

7 

bound  vortex  shroud  distribution 

*lj 

Kronecker  delta 

€ 

shroud  camber 

€ 

e 

effective  shroud  canOjer 

Vi 

X 


4 

P 

o 


T 


T 

A 

T 

0) 


(  )' 


r(r) 

f 


9 


n 


shroud  chord  to  dlsmter  ratio,  c/2R 

propeller  tip  to  shroud  radius  ratio,  R_/R 

? 

c(»stant  fluid  density 
■U/JfVi 

duMiy  tisMi  variable 
Or 

UT/B 

argustent  of  Legendre  function 

denotes  total  differentiation  of  a  function 
with  respect  to  Its  arguiaent 

bound  blade  circulation 

r/RpU 

maxinun  of  f 

nondi&ensional  Fourier  coefficients  of  blade 
circulation 

actuator  disk  lisilt  of  total  aaxinusi  circulation 

difference  between  elesient  and  field  variable; 
e.g.,  Ax  =  {x-\) 

duanv  ang\U.ar  variable  of  integration 

characteristic  functions  giving  effective  shroud 
canber 

constant  angular  velocity  of  propeller 


vii. 


THRBB-DIMENSIONAL  THEORY 
OF 

DUCTED  PROPELLERS 


Since  the  cooc«pclon  of  the  ducted  propeller  or  Kort 
nozzle  in  the  early  1930 coaeiderable  interest  has 
been  aroused  in  this  configuration,  particularly  for  VTOL 
applications  in  recent  years.  Despite  its  proad.se,  however, 
its  inherent  "parsawtriC*  coMplexity  has  precluded  the 
achieveaient  of  coaipzehensive  theoretical  and  experiaental 
developsMnts .  consequently,  the  design  of  an  optimum  ducted 
propeller  has  been  impossible. 

Previous  theories  have  been  reviewed  by  A.  Sacks  and 
J.  Burnell^,  L.  Meyerhof f  and  A.  Flnkelstein^,  and  j.  Megrue^. 
There  are  four  main  categories:  (i)  momentum  methods,  (ii) 
singularity  distribution  methods,  (iii)  "boundary  value"  or 
streasdine  sMthods,  and  ^_v)  combinations  of,  or  approxima¬ 
tions  to,  'chese  methods,  in  general,  though,  the  basic  effects 
associated  with  a  propeller  of  finite  blade  number  are  riaaoved 
by  its  representation  as  a  disk  supporting  a  unifonn  pressure 
jump^'’^'®. 

The  objective  of  the  present  investigation  is  to  derive 
a  three-dimensional  theory  without  the  disk  approximation. 

The  simplest  ducted  propeller  situation  is  choson  %<hlch  retains 
most  of  the  principal  parameters  and  yet  remains  amenable  to 
mathematical  solution. 
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This  case  corresponds  to  the  ducted  propeller  with  discrete 
blading  at  zero  angle  of  attack  In  a  uniform,  unbounded  stream 
of  Invlscld,  Inccmpresslble  fluid.  The  shroud,  blade  emd  hub 
thicknesses  are  assumed  to  be  zero  and  the  shroud  camber  as 
well  as  the  ratio  of  the  blade  chord  to  the  shroud  chord  are 
considered  small. 

The  appropriate  theoretical  model  consists  of  a  cylindri¬ 
cal  surface  oJ  distributed  vortices  positioned  with  Its  axis 
coincident  to  a  uniformly  tr2mslatlng  and  rotating  stream. 

The  blades  are  Introduced  as  Individual,  radial  vortex  lines 
of  varying  circulation  with  accompemylng  helical  vortex  sheets. 
The  blade  circulation  distribution  is  specified  arbitrarily. 

Formulation  of  the  problem  follows  the  general  procedure 
of  thin  lifting  surface  theory.  The  Inherent  periodicity  of 
the  shroud  loading  and  the  harmonic  nature  of  the  oropeller 
field  reduces  the  governing  equations  to  decoupled,  single 
integral  form.  Either  dir<*.rt!  ir  by  modifications,  solutions 
are  expressible  in  terms  of  the  results  of  J.  Weiss Inger^. 
certain  self-sustaining  shroud  load  distributions  appear 
which  seem  2malogou8  to  the  recent  Instability  results  of 
H.  Ludwieg^®. 

It  Is  hoped  that  this  solution  will  provide  not  only  al>so- 
lute  Information  for  design,  optimization  emd  experimental 
test  programming,  but  may  Justify  extended  effort  for  the 
more  difficult  cases  of  angle-of-attack  and  static  thrust. 


CHAPTER  cane 


BASIC  FORMOLATION 

1.1  Phvlol  Probl«»  >nd  Qi>t«1  Background 

lh«  problan  «Ailch  we  have  chosen  to  study  is  s  ducted 
propeller  In  steady  forward  flight  at  zero  incidence  in  an 
incoepressible  fluid  otherwise  at  rest.  The  forward  velo¬ 
city  will  he  assusMd  sufficiently  large  and  the  shroud  casher 
and  blade  loading  low  enough  that  the  theory  of  linearized 
perturbations  uy  be  v^seu.  For  this  investigation  the 
Reynolds  nusiber  will  be  taken  l^jrge  enough  to  confine  vis- 
ccu.«  effects  to  the  boundary  layer,  ibis  will  be  reasonably 
valid  except  in  the  neighborhood  of  the  propeller  tip  if  it 
penetrates  the  shroud  boundary  layer. 

With  respect  to  the  geosMtry  of  the  ducted  propeller 
configuration,  any  thickness  effects  due  to  the  hub,  pro¬ 
peller  blade  or  shroud  will  be  postponed  to  subsequent 
investigations.  However,  the  effective  casi>er  whlcdi  they 
produce  can  readily  be  incorporated  within  the  frawework  of 
the  present  analysis. 

Because  of  the  viscous  skin  drag  associated  with  long 
shrouds,  shrouds  of  Moderate  chord  tc.  diasteter  ratio  are  of 
principal  interest.  Fortunately,  this  is  consistent  with 
the  usual  thin  airfoil  approxlMations  of  satisfying  the 
boundary  conditions  on  the  chordline  by  a  oound  vortex  dis¬ 
tribution  on  the  chordline.  rdchesumn  wd  Weber  have  sug- 
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gestcd  a  value  of  two  as  eui  upper  limit  to  this  ratio  for 
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finite  length  effects  in  their  work  on  annular  airfoils  . 

On  the  other  hemd.  the  ratio  of  the  shroud  chord  to 
the  blade  chord  is  expected  to  be  fairly  large.  That  is, 
for  conventional  propeller  blades,  the  aspect  ratio  is  quite 
large  and  the  shroud  should  be  approximately  of  the  same 
diameter  as  the  propeller.  Therefore  in  order  to  gain  as 
much  thrust  improvement  as  possible  by  preventing  slipstream 
contraction,  the  shroud  has  to  be  lengthened.  Since  this 
increases  the  viscous  shear  drag,  2m  optimum  configuration 
can  exist.  Sufficient  for  the  moment,  it  is  reasonable  to 
take  the  shroud  length  much  greater  than  the  blade  chord  as 
stated. 

In  ordinary  thin  airfoil  theory,  we  have  two  classes 
of  problems:  (i)  the  direct  one  in  which  the  pressure  or 
velocity  on  the  airfoil  is  given  and  the  airfoil  geometry 
is  to  be  determined;  and  (ii)  the  inverse  problem  in  which 
the  airfoil  geometry  is  given  and  the  pressure  or  velocity 
on  the  airfoil  is  to  be  determined.  For  the  case  we  are 
exeuaining  here,  we  have  two  bodies,  the  shroud  and  the  pro¬ 
peller,  and  hence  we  have  four  possibilities  of  combining 
the  above  appropriately,  in  other  words,  we  have  either 
(i)  the  direct-direct;  ^ii)  the  direct-inverse;  >iii)  the 
inverse-direct;  or  (iv)  the  inverse-inverse,  where  the  first 
word  refers  to  the  shroud  and  the  second,  the  propeller.  In 
consideration  of  engineering  applications,  comparison  with 
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experimental  measurements,  configuration  optimization  and 
possible  boundary  layer  studies,  the  Inverse-direct  problem 
has  been  selected.  Nevertheless,  throughout  much  of  the 
analysis,  It  Is  Immaterial  as  to  which  case  is  considered. 

nie  shroud  geometry  will  be  assumed  as  well  as  the  blade 
circulation.  He  are  not  directly  concerned  with  the  chordwlse 
loading  of  the  blade  or  the  determination  of  Its  geometry 
since  we  can  replace  It  with  a  radial  vortex  spike  by  virtue 
of  our  shroud -blade -chord  ratio  assumption.  The  problem,  then, 
will  be  to  obtain  the  shroud  pressure  distribution. 

To  formulate  the  problem,  appropriate  singular  solutions 
of  Euler's  equations  will  be  distributed  to  represent  the 
shroud,  blades  amd  their  shed  vortices.  In  general  the 
strengths  of  these  distributions  will  be  determined  from  the 
requirement  that  the  component  of  the  total  velocity  normal 
to  the  shroud  surface  be  zero  and  that  the  Kutta-Joukowski 
condition  be  satisfied  along  the  shroud  trailing  edge. 

1.2  Coordinate  System 

In  a  space-fixed  coordinate  system,  the  flow  Is  unsteady 
due  to  both  the  motion  of  the  duct  and  the  propeller.  The 
duct  has  the  uniform  translation  U  and  the  propeller,  the 
constant  rotation  H  as  well  as  the  translation  U  ,  see 
Fig.  1.1.  Because  of  this  rotation,  the  motion  is  also 
vmsteady  in  a  uniformly  translating  system  fixed  in  the  duct. 
However,  consider  a  uniformly  tramslating  and  rotating  system 
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fixed  to  the  propeller,  sey  the  Cartetlen  coordinetes  (x,y,z} 
or  the  equivalent  cylindrical  coordinates  (x,r.6)  .  Increas¬ 
ing  B  is  chosen  in  accordance  with  the  righthand  rule  for 
positive  advance  along  the  positive  x-axis  and  is  sMasured 
fron  the  positive  y-axis. 

In  this  systSM,  Fig.  1.2,  we  see  that  the  inconing  "free 
streaa"  has  an  axial  coaqponent  of  Magnitude  U  in  tho  posi¬ 
tive  x-direction  and  a  tangential  cos^onent  flr  in  the  posi¬ 
tive  or  increasing  d -direction.  How,  the  blades  are  at  rest 
With  respect  to  the  propeller  coordinates,  but  the  duct  is 
rotating  with  angular  speed  o  about  the  propeller  axis  as 
shown.  Since  the  duct  is  axisysaetrical  idMut  its  axis  of 
rotation  and  only  the  nornal  coaponent  of  the  velocity  is 
required  to  vanish,  the  rotational  notion  of  the  duct  does 
not  affect  the  potential  solution.  In  other  words  this 
action  is  arbitrary  and  can  be  set  equal  to  zero  in  the  pro¬ 
peller-fixed  coordinates.  This  nathanatically  equivalent 
Model  t^xen  corresponds  to  the  physical  case  in  «diich  the 
shroud  is  rotating  along  with  the  propeller,  see  Fig.  1.3* 

A  sisgile  case  night  Illustrate  this  nathenatical 
equivalence,  consider  a  tvo-dinensional  thin  airfoil  in  a 
unlfom,  inviscid  strean.  Then,  suppose  an  infinitesiaally 
thin  belt  is  stretched  around  the  airfoil  emd  set  into  Mtion. 
For  the  resultant  flow,  the  nomal  boundary  and  Xutta  con¬ 
ditions  arc  unaltered  by  this  notion  and  so  the  solrtion  is 
the  sane  as  before,  in  a  real  fluid,  of  course,  the  boundary 
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layer  would  be  changed  and  hence  also  the  potential  flow, 
proportional  to  the  modification  of  the  effective  airfoil 
shape.  This  correction  Is  of  higher  order,  however,  and  so, 
consistent  with  the  usual  thin  airfoil  approximations,  the 
two  are  equivalent. 

If  the  axis  of  the  duct  were  Inclined  at  an  angle  of 
attack  u  with  respect  to  the  translational  motion,  no 
equivalent  steady  flow  would  exist.  That  Is,  In  this  case 
there  are  tlaw  dependent  Incoming  velocity  components.  With 
the  stream  U  coincident  with  the  xs-plane  at  tine  t  «  0  , 
we  have  for  the  free  steam  velocity  3^  In  propeller-fixed 
coordinates , 

q^  *  U  cos  a  +  ftr 

+  U  {  sin  (a-Ot)  ij.  + 

cos  (©-fit)  1^  ]  sin  a  (1«1) 

setting  a  3  0  ,  we  see  Eq.  (1-1)  reduces  to  the  steady  com¬ 
ponents  u  and  Hr  we  have  already  noted. 

1.3  The  Mathematical  Model 

with  these  considerations  w«  can  proceed  to  a  de¬ 
tailed  examination  of  the  model.  Fig.  l.^(  gives  the 
geometrical  properties  of  the  configuration.  The  plane 
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determined  by  the  leading  edge  of  the  duct  is  perpendicular 

to  the  x-axls  at  x  =  -c/2  auid  the  trailing  edge,  at 

X  =  ►c/2  .  The  propeller  disk  of  radius  is  normal  to 

the  x-eocls  at  x  =  x^  with  the  of  the  N  blades 

P 

lying  in  the  xy-plane.  The  slope  of  the  shroud  with  respect 
to  the  x-axis  in  any  xr -plane  is  designated  by  tan  €(x)  , 
which  is  taken  positive  for  increasing  shroud  radius  with 
Increasing  x  .  A  reference  shroud  radius  is  R  .  For 
convenience  with  respect  to  tip  clearance,  R  is  taken  as 
the  radius  of  the  shroud  caniber  surface  in  the  propeller  plane. 

In  accordamce  with  our  assumptions,  we  can  now  repre¬ 
sent  the  configuration  by  a  distribution  of  bound  auid  shed 
vortices.  First,  the  propeller  model  follows  classical 
theory,  nie  circulation  distribution  of  the  bound  radial 
vortex  spikes  replacing  the  blades  is  designated  by  r(r^)  , 
the  subscript  v  being  generally  used  to  designate  the 
location  of  a  vortex  element.  Each  blade,  by  symmetry,  has 
identical  loading  and  is  accompanied  by  its  shed  semi-infinite 
helical  vortex  sheet  of  strength  df/dr^  per  unit  radial 
length.  The  helical  shape  is  assumed  to  be  determined  by 
the  Incoming  stream  components  at  r^  ,  its  point  of  shedding. 
For  the  shroud  we  have  a  form  of  a  lifting  surface.  The  bo\ind 
vortex  distribution  per  unit  length  is  >{x^,6^;  r^)  where 
r„  =  r  (x„)  denotes  the  radius  of  the  shroud  caniber  surface. 
Because  of  -y's  angular  dependence  generated  by  the  finite 
blade  number,  each  vortex  ring  has  a  trailing  vortex  cylinder 
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made  up  of  an  infinite  nuniber  of  semi-infinite  helical  vor¬ 
tices  of  strength  ( )  ( r^dS )  .  The  positive  senses  of 
these  vortices  will  be  defined  later. 

For  a  given  propeller  loading,  we  can  now  formulate  an 
equation  to  determine  y  .  Let  3  be  the  total  velocity 
vector  given  by 


S  =  %  +  3r  +  Sr'  +  Sy  +  Sy' 


(1.2) 


where  3^  has  been  defined  previously  and  3p  is  the  velo¬ 
city  induced  by  the  N  boxmd  blade  vortices;  ,  by  the 

shed  blade  vortex  sheets;  ^  ,  by  the  shroud  hound  vortices; 
and  Sy '  '  shed  shroud  vortices.  With  3  •  i  and 

S  ‘  ibj.  A*  the  axial  and  radial  components  of  q  respectively, 
then  for  the  flow  to  be  everywhere  tamgent  to  the  shroud,  we 
require  from  Fig.  l.t. 


tan  €(Xg) 


(1.3) 


vrherc  r_(x_)  has  been  defined.  Eq.  (1.3)  is  exact  within 
our  model  limitations.  If  we  assume  e  «  7r/2  ,  tarn  €  =  e+0(€ 
'Ihus  for  a  moderately  loaded  propeller  ox  0(r/RpU)  «  1  , 

3  •  jL  =  0(u)  and  hence  from  Eq.  (1.3).  3  *  ij-  “  ®("^'  ”  0(U*;) 

This  leads  to  the  linearization  of  the  RHS  of  Eq.  (I.3).  Expan^ 
ing  q  •  and  3  *  A  as 
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^•i 


"ST 


(r,-S)  +  ... 


r,=R 


5q  •  ^ 
+  -^?— 


+  •  • 


r^=R 


(1.^) 


amd  noting  that  (r^'R)  =  0(c€)  ,  we  have  to  retain  only  the 
first  terms  from  consideration  of  the  previous  orders.  In 
addition,  from  Eq.  (1.2)  with  a  =  C  for  ^  of  Eq.  (l.l), 
the  axial  component  can  be  replaced  by  U  and  Eq.  (1.3) 
reduces  to  a  familiar  thin  airfoil  type  approximation, 


3(x,.R.eg)  •  4 
u 


(1.5) 


That  is,  we  have  only  to  satisfy  a  perturbation  boundary  con¬ 
dition  on  a  mean  shroud  or  reference  surface.  Since  the  radial 
component  of  the  free  stream  is  zero  for  zero  angle  of  attack, 
Eq.  (l.l),  we  have  only  to  calculate  the  Induced  radial  com¬ 
ponents.  The  bound  emd  sbed  shroud  vortices  as  formulated  are 
located  by  r^  .  But  consistent  with  the  orders  we  have 
retained,  an  expansion  as  carried  out  for  £  •  i,  and  q  • 
shows  that  use  of  the  meam  radius  R  is  sufficient.  In  sub¬ 
sequent  discussions,  q^  and  ^  •  will  be  so  considered 
without  further  notation. 
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Eq.  (1.5).  together  with  a  Kutta  condition  at  the  trail¬ 
ing  edge  of  the  shroud,  is  the  basic  equation  for  the  deter¬ 
mination  of  the  shroud  vorticity  and  hence  the  loading.  Hie 
propeller  inflow  could  also  be  calculated  and  hence  the 
required  blade  geometry  for  the  selected  blade  loading. 

1.4  Velocity  Influence  Functions 

In  order  to  express  the  radially  induced  flow,  or  wash, 
of  Eq.  (1.3)  as  well  as  the  inflow  components,  it  is  conveni¬ 
ent  to  introduce  three  elementary  velocity  fields.  These 
fields,  or  influence  functions, 
of  unit  strength  and  tult  length  which  lie  in  the  axial,  radial 
and  circumferential  directions  respectively  as  shown  in  Fig.  1.5 
In  general,  then,  we  have  from  the  Biot-Savart  law  for 
the  velocity  induced  by  an  infinitesimal  vortex  in  an  incom¬ 
pressible  fluid, 

1  I  ^  D 

Jh-  “  I?  fjp 

where  I  is  a  unit  vector  in  the  positive  direction  of  our 
desired  element  and  O  is  the  vector  from  the  element  at 
^*v'^v'®v^  to  the  field  point  (x,r,9)  , 


arc  due  to  vortices 


D  =  (x-x^)  i  +  (y-y^)  I  +  (z-z^)  H 


(1.7) 
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with  the  relatiorj 

X  “  X 

y  =  r  cos  6 

*  =■  r  sin  0  (1-8) 

between  the  cylindrical  and  Cartesian  coordinate  systems.  The 

unit  vectors  ^md  k  are  directed  along  the  x,y  and  z-eoces 

respectively. 

From  Fig.  1.5  we  have  for  the  axial  element,  or  , 

simply 

1=1  (1.9) 

where  positive  circulation  is  chosen  in  accord  nee  with  a 
righthand  screw  advancing  in  the  positive  x-direction.  For 
the  radial  element,  or  ,  again  using  the  righthand  rule, 
we  have 

I  =  cos  2  ®v  —  (I.IC) 

with  positive  circulation  corresponding  to  advauice  in  the 
positive  radial  direction.  Finally,  for  the  circumferential 
element,  or 

=  -sin  0^  2  +  cos  0^  k 


I 


(1.11) 
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for  positive  circulation  in  the  increasing  0-direction.  This 
sign  convention  on  the  shroud  will  turn  out  opposite  to  that 
usually  esiployed  in  wing  theory. 


1.5  Calculation  of  the  Induced  Velocities 

From  Jlp  and  ^  ,  the  velocity  fields  for  both  the 
propeller  and  shroud  vortices  may  be  simply  expressed, 
and  from  ^  and  ^  ,  the  fields  of  the  shed  helical 
vortices. 

For  .  we  Integrate  over  each  blakde  and  sum  the 
results  over  all  blades.  Or,  from  Fig.  1.4,  we  have 


z 


1=1  c 


(1.12) 


where  Eqs.  (1.6)  through  (1.8)  together  with  Eg.  (1.10)  give 

1  2ind  the  argument  of  I  locates  a  vortex  element,  the 
~P  ~P 

field  point  remaining  tinspecifled  for  the  nxment.  the  index 
t  denotes  the  blade  nxiraber,  the  blade  being  chosen 

arbitrarily  to  coincide  with  the  y-axls. 

Likewise  for  the  shroud,  we  have  by  integration  over  y 
placed  on  the  reference  cylinder  of  radius  R  , 


,>c/S 


^  is  given  by  Eqs.  (1.6)  through  (1.8)  and  Eq.  (1.11). 
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To  express  qj,^  and  3^'  *  require  the  velocity 
Induced  by  a  semi -infinite  helical  vortex  of  unit  strength, 
with  its  eocis  coincident  with  the  x-axis.  As  assumed,  the 
shape  is  determined  by  the  free  stzeam.  That  is.  the  coor¬ 
dinates  of  2my  point  on  the  helix  which  originates  at 
(Xy.r^.0^)  are  given  by  ( x^+Ut , r^ , )  .  Here  x  is 
the  conventional  variable  representing  the  time  for  a  vortex 
element,  convected  by  the  free  stream  to  travel  from  the  point 
of  shedding  to  that  point.  Subsequently,  then,  we  recognize 
that  the  helical  filament  can  be  represented  by  an  Infinite 
series  of  component  elements  parallel  and  concentric  to  the 
x-axls.  Integrating  over  these  component  elements  of  length 
Udx  and  r^^ldx  respectively,  we  obtain  the  desired  velocity 
field.  This  can  be  compared  with  the  results  of  other  pro¬ 
peller  investigations,  see  T.  Moriya^^,  for  exanqple. 

From  Helmholtz's  theorem,  the  strength  per  unit  length 
of  the  vortices  shed  from  the  propeller  is  -P'.  The  prime 
denotes  differentiation  with  respect  to  the  argument  and  the 
adopted  sign  convention  introduces  the  minus  sign.  Carrying 
the  integration  and  summation  over  the  propeller,  we  have 

^  0 

Where 


=  Ijj(Xp+Ux,r^.27r//Srnx) 
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=  I^(x^+UT,r^.27ri/N+nT)  (1.1^) 

is  given  by  Eqs.  (1.6)-(1.8)  auid  Eq.  (l.9)i  emd  was 

given  for  ^  .  The  integration  over  the  dummy  time  parameter 
T  is  identified  as  the  velocity  field  of  the  semi-infinite 
helix  of  unit  strength. 

Finally,  we  cam  derive  ^ »  for  the  shroud,  where  the 
strength  per  unit  length  of  the  shed  vortices  is  . 

Integrating  over  the  reference  surface  of  radius  B  as  for 

%  ' 

*>'  't/2  L 

where, 


I,  2  ^,(x +uT,R,e +n'r)  (1.15) 

“U  ~~C.  V  V 


and  were  previously  noted. 

Superposition  of  the  velocity  fields  expressed  by 
Eqs.  (1,12),  (1.13).  (1-1^)  and  (1.15)  determine,  within  the 
limits  of  our  approximation,  the  total  induced  velocity  at  aaiy 
field  point. 


and  1^ 
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1.6  Fonrolatlon  of  B«»lc  Equation  for  y 

With  the  total  induced  velocity  field  we  can  complete 

the  explicit  formulation  of  Eq.  (l.i>)  for  the  shroud  vortex 

distribution  y  .  in  particular,  we  need  only  the  radial 

conponents  of  the  velocity  evaluated  on  the  reference  surface. 

TO  determine  these  components,  we  first  carry  out  the 

required  scalar  triple  products  •  1.  D  .  ihe  radial  unit 

vector  i^  is  given  by  Eq.  (1.10)  with  0^  replaced  by  9  r 

i  ,  by  Eqs.  (1.9),  (I.IC)  and  (l.ll);  and  D  ,  by  Eqs,  (1.7) 

and  (1.8).  With  x  =  x  ,  r  =  R  and  6=9  ,  we  have 

8  8 


^  X  D  =  -  r^  Sin 

ir  •  Ip  ^  D  »  -  (x,“Xv)  ‘in  (0,-6^) 

i^  •  X  D  =  (x^^x^)  cos  (0s-®v) 


(1.16) 

(1.17) 

(1.18) 


The  other  quemtity  we  need  is  the  magnitude  of  the  vector  D 
from  any  element  to  the  field  point  given  by 

=  Ax^  +  +  r^  -  2Rr^cos  A0^  (1.19) 

for  field  points  on  the  shroud  where  Ax^  =  (x^-x^)  and 
A0  =  (0  -0  )  . 

The  radial  shroud  velocity  components  now  follow  from 
Eqs.  (1.12),  (1.13),  (1.1^)  emd  (1.15).  For  the  propeller 


contributions,  the  bound  vortices  give. 


Sr 


i_  * 


i=i 


■  rrZ,® 


r  p  _3 

in  D  ^  (1-20) 


where  Ax^  h  (x^-x^)  ,  A0£  =  (0^~2Tr£/f!)  and  D  is  given  by 
Eq.  (1.19)  with  x^  replaced  by  x^  and  0^  ,  by  2n£/N  . 
For  the  shed  vortices. 


1  r”p  r 

•  -r  *  J  V  j  “i" 


Ae 


i-l  0 


iT 


-  HAXp^  cos  AO^^)  d"^  drdr^  (1.21) 

with  AXp^  H  '  ^®iT  =  (ea-2iTi/N-flT)  ;  and 

replaced  by  (Xp+ur)  and  0^  by  (gxi/N+Ox)  in  D  .  For 
the  bound  shroud  vortices,  we  have 


5y 


c/2 

7t 

4  “ 

1 

-c/2  ^ 

f  VD-3 

cos  A0^ 

replaced  by  R 

in  Eq. 

(1.19);  and 

(1.22) 


--C/2  tt  oo 

3»-  •  ix  ■  "VT 


-  CO.  09^^)  d‘3  dx  dS^d*^ 


(1.23) 
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for  Che  shed  shroud  vortices  with  ' 

^®vT  ^  \  replaced  by  (x^+Ut)  ,  r^  by  R 

and  6^  by  in  Eq.  (1.19). 

If  Eqs.  (1.20),  (1.21),  (1.22)  and  (1.23)  together  with 
Eq.  (1.19)  are  substituted  into  Eqs.  (1.2)  imd  (1.5)*  the 
basic  equation  for  the  shroud  vortlclty  is  determined  in  terms 
of  a  surface  or  double  integral  relation.  As  such  it  has  the 
general  form  familiar  to  lifting  surface  theory.  That  is,  for 
a  specified  propeller  loading,  the  propeller  terms  may  be 
grouped  with  the  shroud  camber  term  forming  ui  effective  camber. 
The  shroud  bound  and  shed  vortices  are  then  required  to  produce 
a  flow  tangent  everywhere  to  this  fictitious  surface,  as  well 
as  to  satisfy  a  Kutta-Joukowski  condition  at  the  trailing  edge. 

Not  only  does  this  equation  have  this  general  form,  but  it 
may  be  compared  in  detail  with  the  ring  wing  surface  theory  of 

Q  1  p 

J.  Weissinger  .  In  fact  if  fl  is  set  equal  to  zero  and  the 
T  integration  carried  out,  the  form  is  identical  with  that  of 
a  modified  ring  wing  having  non-axlsymraetrlc  camber. 

On  the  other  hand,  if  the  blade  loading  is  not  given,  then 
a  second  equation  must  be  written.  In  view  of  our  approxima¬ 
tion,  a  "lifting”  line  iu>del  would  be  appropriate.  This  equa¬ 
tion  would  be  coupled  with  Eq.  (1.5)  by  the  ch2mge  in  the 
effective  local  angle  of  attack  of  each  blade  element  pro¬ 
duced  by  the  total  inflow. 
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DERIVATION  OF  GOVERNING  EQUATIONS 

2.1  Periodic  Expansion  of  Shroud  Vortex  Distribution 

In  general  the  inversion  of  a  lifting  surface  type  equa¬ 
tion  Is  exceedingly  complex.  However,  in  the  present  case  as 
well  as  for  a  ring  wing,  the  cloning  of  the  shroud  surface 
introduces  a  natural  periodicity  which  yields  an  essential 
simplification. 

In  particular  the  str-  ngth  of  the  bound  shroud  vortices 
may  be  expanded  in  a  complex  Fc  ' i er  series  to  permit  inte¬ 
gration  over  the  0^  variable,  ct 

00 

>(x,0)  H  u  ^  cjx)  e^®  (2.1) 

m=-® 

The  speed  u  is  inserted  for  conven  ,ence  to  non-dimsnsion - 
alize  certain  subsequent  expresslonr ,  Only  mN  components 
are  used  as  indicated  by  the  periodicity  in  27r/N  of  the 
propeller  velocities  of  Eq.  (1.20)  and  (i.2l)  for  the  equal 
blade  loading  of  the  present  case.  Since  y  is  real,  we 
require  that  the  complex  coefficients  for  positive  (+m) 
and  negative  (-m)  integers  be  complex  conjugates,  or 


2C  „  ^  +  i  A„ 

-m  m  m 


^^+m  *  ®m  ~ 


(2.2) 
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where  A  and  B  are  the  real  coefficients  of  the  sine  and 
m  m 

cosine  terms  respectively  of  a  trigonometric  Fourier  series  over 
positive  m  for  7  . 

For  m  =  0  ,  we  note  Aq  =  C  and  2Cq  =  Bq  ;  i.e., 

UBq/2  is  the  zeroth  Fourier  component  of  the  shroud  vorti- 
city  in  the  propeller-fixed  coordinates,  or  equivalently,  the 
time-independent  part  in  the  shroud-fixed  coordinates.  As 
such  its  physical  role  will  be  shown  to  determine  the  time 
average  difference  in  shroud  inner  and  outer  static  pressures. 

2.2  Reduction  of  Shroud  contributions 

Before  substitution  of  Eq.  (2.1)  into  Eqs.  (1.22)  and 
(1.23)  for  the  shroud  contributions  to  the  radial  velocities, 
we  introduce  the  radius  R  as  a  reference  length,  or 

X  2  x/R  ,  . . .  (2.3) 

This  will  limit  the  explicit  occurrence  of  the  shroud  chord 
to  diametei  ratio  defined  as 

(S.O) 

to  the  limits  of  integration  over  the  coefficients  Aj„  and 
.  We  also  change  our  variable  of  integration  from  9^  t.o 
Ae^/2  .  Since  the  integrand  is  periodic  in  0^  ,  we  adjust 
the  new  limits  of  integration  from  -tx/2  to  ■n/2  ■  In  add;.- 


tion  ,  the  separation  distance  D  is  simplified  by  means  of 
the  cosine  double  angle  formula. 

Now  using  Eq.  (2.1)  in  terms  of  the  new  quantities  in 
2q.  (1.22),  the  order  of  integrations  and  svunmation  may  be 
interchanged  because  has  at  worst  a  square  root  singu¬ 

larity  at  the  shroud  leading  edge.  We  can  then  split  out 
the  0^  dependence  from  the  integrations  to  obtain, 


3,  •  ir  “  “I 

m=  -a> 


(2.5) 


where, 


■•ym 


Ax. 

~2 


-  [ 


V2 


-2imN9 


cos  28 


‘-n/2 


(Ax^^  t  4  sin^ 


d0 


(2.6) 


In  other  words,  Eq.  (2-5)  represents  the  Fourier  expansion 

of  radial  velocity  induced  on  the  shroud  by  the  bcund  shroud 

vortices.  The  coefficients  of  this  expansion  are  si^v-^-Y  the 

integrals  of  the  coefficients  of  the  shroud  vortices  operated 

on  by  the  functions  of  Eq.  (2.6). 

These  functions  may  be  simplified  and  expressed  in  terms 

o2  known  functions.  The  sine  term  of  the  exponential  is  odd 

2 

and  cos  28  and  sin  8  ,  even,  so  that  the  imaginary  part 
of  is  zero.  For  the  real  part  we  first  note  the  trig¬ 

onometric  identity 
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2  cos  2inN0  .  cos  29  = 

C03  2(niN+l)0  +  cos  2(mN-l)0 


(2.7) 


and  then  use  the  trigonometric  integral  form  of  the  Legendre 

fxinction  of  the  second  kind  and  half  order  Q  .  ,  developed 

n-sj 

l-a 

under  this  program 


-  J 


jr/2 


cos  2n9 


-ii/2  (2(a>-l)  +  k  sin  9) 


d0 


(2.8) 


by  differentiation  with  respect  to  the  parzuneter  ®  .  Com¬ 
paring  the  resulting  forms,  we  have 


“Vm 


^v 


mN 


(\) 


s  (1  +  Ai^Va) 


(2.9) 


^mN  -  ^Wl/2  *  °mN-3/2)  denotes 

total  differentiation. 

Alternatively,  may  be  expressed  in  complete  ellip¬ 

tic  integrals  or  Riegels*  functions^^.  For  the  higher 
harmonics,  the  Legendre  functions  are  more  convenient.  For 
any  order  they  have  the  same  logarithmic  singularity  at  an 
argument  of  unity  and  decrease  raonotonically  to  zero  for 
increasing  argument,  see  Fig.  2.1.  At  any  argument,  they 

1  ii 

decrease  monotonically  with  order,  NBS  tables  ,  as  well 
15 

as  extended  tables  prepared  'inder  this  study  are  available. 
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TO  derive  an  equation  for  the  vortices  shed  from  the 
ehroud  in  a  form  like  Eq.  (2-5),  we  can  proceed  along  similar 
lines.  The  only  essential  complication  appears  in  the  addi¬ 
tional  integration  over  the  dummy  time  variable  t  .  There¬ 
fore,  we  Interchange  the  order  of  integration  of  t  wd  9^ 
and  then  use  dd^^/2  as  the  variable  of  integration  instead 
of  6^^  .  This  gives  a  corresponding  in  terms  of  an 

integration  over  the  Legendre  functions.  That  is,  we  find 

^  I  ^  m(^v^  ^  (2. 1C) 

m«-«»  -X 

where  P  '  is  given  by  a  Fourier  integral  over  the  dimen- 

#  iO 

sionless  dummy  tine  variable  t  =  fn  ,  or 


p  ,  s  jyss.  r  ^“iwN^  j  ^  g'  (3)_) 
S  m  TwT  ^  “^VT  laN'  2^ 


ajg  =  1  + 


(2.11) 


^nN  =  ^Wl/2  •  °m-3/2^  •  propeller 

advzmce  ratio  J  and  the  propeller  tip  clearance  parameter 


U  defined  by 


?8 


J  ^  U/fiRp  (2.12) 

It  =  »lp/R  (2.13) 


represent  the  "shroud  advance  ratio"  U/flR  .  The  factor 
imN  results  from  the  differentiation  of  the  shroud  vor- 
ticity.  This  differentiation  which  is  carried  out  term- 
wise  is  justified  on  a  physical  basis. 

2 . 3  Decomposition  of  Governing  Equation  for  7 

The  introduction  of  the  expansion  of  Eq.  (2.1)  for  y 

has  permitted  the  explicit  integration  over  6^  .  if  we  now 

express  the  propeller  contributions  in  terms  of  their  Fourier 

components,  we  can  decompose  the  single  equation  for  the 

vortex  distribution  over  the  shroud  surface  into  an  infinite 

set  of  uncoupled  equations  integrated  along  the  chord. 

Let  UCr  (x  )  and  UCp'  (.x  )  be  the  complex  Fourier 
1  m  s  I  m  s 

coefficients  of  the  radial  velocity  Induced  on  the  shroud  by 
the  bound  and  shed  propeller  vortices,  or 


in=  -«> 


( 2  - 14  ) 


qp'  ' 


1 


u 


(2.15) 
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To  find  these  coefficients  we  have,  from  orthogonality,  the 
relations 


Tm 


1 

2nv 


/-TT 

1 


^-iinNe, 

e  s 


d^ 


(2.16) 


1  d% 


(2.17) 


for  a  period  of  2’’^  .  Since  •  i^  and  fip'  ’  A,,  real 

quantities,  the  values  of  these  coefficients  corresponding  ro 
+m  end  -m  ,  must  be  complex  conjugates . 

First,  substituting  Eq.  (1.2C)  into  Eq.  {?.l6),  we  trans¬ 
form  to  the  non-dimensional  quantities  previously  used  and 
interchange  the  order  of  the  r^  ”  and  6g  ”  integrations. 

SiP'-e  the  0^  dependence  is  periodic  over  the  range  of  inte¬ 
gration,  we  adjust  the  limits  to  (-tt  ^  ?7r£/N)  ^  ^  (^  ‘‘^  2r£/ 

and  integrate  on  d  =  Ad^/2  for  -7t/2  ^  9  ^  f^/2  .  The  depend¬ 
ence  over  the  blade  summation  index  £  ,  then,  disappears  from 
this  integral  since  it  occurs  only  in  the  exponential  power, 
-2''nn^i  ,  which  is  equal  to  one  for  any  integer  mi  •  Tliat  is, 
for  N  identical  blades,  each  blade  contributes  equals/  to 
any  harmonic.  The  integral  on  9  itself,  after  removing  a 
factor  of  from  5^  ,  is  evaluated  in  terms  ^'f  the 

Legendre  functions  ^7  means  of  an  equation  analogous 

to  Eq.  (2.7)  for  th''  sine  function,  as  v/ell  as  Eq.  (2.3). 
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The  final  result  is. 


]/2i^  (2.18) 


The  non-diiuensional  blade  circulation  r  is  defined  by. 


r  H  r/RpU  (2.19) 

Instead  of  u  ,  Ar  could  be  used,  but  the  translational 
P 

i,  ad  is  better  for  prinary  emphasis  on  the  shroud. 

If  we  substitute  Eq.  (1.21)  into  Eq.  (2.17),  can  pro¬ 
ceed  in  the  Same  way.  'nie  complication  of  the  integration  over 
the  dummy  time  vsuriable  T  can  be  handled  by  an  integration 
interchange  corresponding  to  the  analysis  for  the  vortices  shed 
from  the  shroud.  Carrying  this  out,  we  find  that 


87t‘ 


r  f'(r-^)  r 


e-^^  {  dip,  S'^(S,) 


-  iJH 


5^  S  1  +  (  diip/  +  J/2i^ 


(2.20) 


Which  is  a  Fourier  integral. 
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In  a  broader  sense,  Eqs.  (2.18)  and  (2.2C)  represent  the 

harmonics  of  the  radial  velocity  of  a  free  propeller  in  an 

unbounded  imiform  stream.  That  is,  ,  R  are  interpreted 

s 

as  the  coordinates  of  an  arbitrary  field  point.  S.  Tsakonas 
and  J.  Breslin  have  determined  the  harmonic  components'^  of  the 
pressure  field  for  a  free  propeller  using  the  same  vortex 
representation.  Their  results  are  expressed  in  complete 
elliptic  integrals  instead  of  Legendre  functions. 

Now,  we  substitute  Eqs.  (2.5).  (2.1C),  (2.l4),  (2.15). 
(2.18)  and  (2.20)  into  Eq.  (1.5)  a^id  equate  the  coefficients 
of  each  harmonic.  After  transposing  the  propeller  terms,  we 
obtain 

*0m  "(^s)  -  Cfm^^p)  -  m('^p)  = 

X 

where  is  the  Kronecker  delta  and  m  extends  over  all 

negative  and  positive  i»itcgers.  The  shroud  camber  appears 
only  in  the  zeroth  harmonic  since  it  is  axisymmetric  by 
assumption.  In  other  v/ords,  in  inviscid  flow; 

The  effect  of  axisymmetric  shroud  camber  for  a  ducted 
propeller  at  zero  angle  of  attack  appears  only  in  the  steady 
shroud  load  for  a  fixed  blade  loading. 
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2.4  Splitting  of  Governing  Complex  Equation 

In  complex  form,  Eq.  (2.21)  constitutes  a  doubly  infinite 
set  of  equations  over  •all  m  for  the  coefficients 
which  determine  the  shroud  vorticity.  Alternatively,  we  can 
separate  the  real  and  imaginary  terms,  producing  an  infinite 
set  of  coupled  pairs  of  equations  in  and  over 

m  =  1,2,3  •• • 

For  simplification,  we  note  from  the  recursion  relation¬ 
ship  for  the  Legendre  functions  , 


(2.22) 


that. 


®mN  “  ^-mN 

T  n:  -T 

niN  -raN 


(2.23) 


From  Eqs.  (2.9)  and  (2.11),  then,  we  have  that 
F^  are  complex  conjugates  and  Fy  ^  and  F^'  are 
complex  conjugates.  That  is. 


7  .-ro 


F  =  v<(F  ) 

7,+m  '  7m' 


F  =  *(F,  ) 
7,-m  '  7m' 


V,+n>  = 


^Vm  =  m)  - 


(2.2i  ) 
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In  addition,  the  propeller  Fourier  coefficients  for  +in  and 

-m  ue  iiompXex  conjugates,  in  fact,  if  Eqs.  (2.23)  used 

in  Eqs.  (2.l8)  and  (2.20;,  tho  equations  for  Cj.  and 

Cp  ^  are  the  sane  as  for  P^,  and  F  _  with  the  real 
* , ~n  y , +n  y , -n 

part  replaced  by  tb'j  inag inary  part,  and  the  equations  for 
Cp'  and  Cp'  are  exactly  the  sane. 

these  results  and  Eqs.  (2.2),  we  can  add  and  sub¬ 
tract  ..  (2.21)  into  a  set  of  coupled  equations 

for  2md  . 

2  fen  '  -  2«(=r'.)  ’ 

‘  \  (2.S5) 


and. 


2ji(Cpjj|)  -  25(Cp'jj)  = 

-X 

where  n  »  0,1,2  ...  .  For  m  --  0  ,  =  0  as  noted  before 

and  Eq.  (2.25)  determines  Bq  . 

The  Intermediate  form  of  Eqs.  (2.25)  and  (2.26)  nay  be 

13 

further  simplified.  First,  from  the  recursion  relation  , 


X 
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^+*5  ”  ®n-3/2 


(2.2,0 


the  sin  inNT  term  of  ,  corresponding  to  the  elements 


convected  by  the  axial  flow,  may  be  reduced  to  Q^_jj  •  On  the 
other  hand,  the  term  produced  by  the  swirl  depends  on  Ax  , 


but  this  is  a  perfect  differential  in  t  ;  i.e.. 


(2.28) 


Integrating  this  contribution  to  by  parts  imd  not¬ 


ing  Q  ,  goes  to  zero  at  infinity,  we  get  a  contribution 
n-’s 


independent  of  t  from  the  lower  limit  and  cjui  add  the  remain¬ 


ing  integral  to  the  reduced  axial  terms. 


The  term  S  (Cp'^^j)  yields  to  analogous  steps  while  X{Cj^) 
requires  only  Eq.  (2.27).  If  we  set  ^-(SCp'^ji)  =  « 


<i-(2Cp„)  =  -  Aj,  and  i(2Cp.  )  =  -  Ap .  ,  we  get 


m  =  J2j  ^  ®mN(“3)  " 


g  „(Ax  ,r  )  )  dr 
^mN'  p  v'  V 


(2.29) 


umN'^Ax-  r  d 


''r-„  ■  f'  i 


(2.31) 


The  functions  and  are  Fourier  integrals  given  by: 

®  "/o  1 

00  2 

■W  “  1  Wl,  ♦  ]  «•  (2-32) 

with  Qnjjj.jj  3md  having  the  argument  3^  . 

Turning  to  the  kernels  of  Eqo.  (2.25)  and  (2.26) »  we  find 
that  they  can  be  expressed  in  simil2ur,  but  simpler,  terms.  In 
fact,  we  can  show  that  one  is  proportional  to  the  derivative 
of  the  other,  or, 

(2.33) 

Therefore,  we  can  define  a  single  kernel,  say  =  (J(i/mN)  • 

•  and  obtain 

V^v)  =  ^  ^  W"l)  "  ^mN^^v'^)  ^  (2.34) 

where  the  argvunent  of  and  of  g^^(Ax^,l)  is 

Sg  because  r^  is  unity.  '^^5^(3^^  )/47r  is  readily  identified 
as  Si  (f,  )  or  the  contribution  of  the  bound  vortices. 

From  Eqs,  (2.29)  through  (2.34),  the  intermediate  forms 
of  Eqs.  (2.25)  and  (2.26)  assume  the  final  forms. 
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2^0in  ^  ■  ®r  m  ■ 

f  -  J  (2-35) 

*“X 

^Fm  '*'  ^r"ra 

f  /  »«'s.‘*^v  ^  /^  's.’s;'*^.  (2-36) 

•“X  “*X 

where  e  is  a  function  of  ;  Ap^  ,  Ap  ,  and  Bj,  ,  of 

\  \  ‘  \  '  \ 
m  =  0  ,  =  0  and  Eq.  (2.35)  determines  Bq  . 

•6  Discussion  of  coupled  Equations 

Eqs.  (2.35)  and  (2.36)  for  the  bound  chordwise  di' 
tions  of  the  shroud  vorticity  contain  several  importcuit  f6ai.wjres 
First,  the  propeller  contributions  depend  on]y  on  the 
axial  separation  of  the  shroud  point  and  the  propeller  plane 
for  a  given  number  of  blades,  blade  loading  and  tip  clearance. 
Since  they  do  not  depend  on  the  eujgular  position  of  the  blade, 
the  effect  of  finite  blade  chord  could  be  incorporated  In  the 
present  formulation.  That  is,  the  blade  could  be  represented 
by  a  distribution  of  radial  elements  with  6^  -  ©^(x^)  eind 
r  =  .  The  corresponding  propeller  Fourier  coeffi¬ 

cients  would  then  be  integrated  over  the  appropriate  range 
of  Xp  . 
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Second,  the  coupling  of  the  equations  is  due  to  the  swirl, 
or  n  ,  of  the  Incoming  flow.  In  other  words,  from  the  inter¬ 
mediate  forms  of  Eqs.  (2.25)  and  (2.26),  the  equations  would 
be  Independent  only  if  were  zero.  But  from  Eq.  (2.11 

this  requires  that  fl  0  .  To  see  this,  we  resubstitute 
T  =  fit  and  use  the  alternative  dimensionless  time  t  =  Ux/R  . 
Then  the  exponential  factor  beccxnes  one  as  fl  -»  0  and  the 
swirl  terms  disappear,  giving  the  real  quantity 


£  im 

n  -*  0 


_  mu 

'•y 'm  ~  4iv 


0 


..  ,  -  .  2 

where  =  1  +  (Ax^-x)  /2  .  The  decoupled  equations  have  the 
same  kernel  which  has  a  Cauchy  singularity. 

Third,  various  singularities  will  appear  in  our  governing 
equations.  In  general  the  propeller  contributions  will  remain 
finite  since  u  <  1  .  For  the  limiting  case  of  p,  -►  1  ,  a 
special  study  of  the  anticipated  behavior  of  the  bound  blade 
circulation  near  the  tip  is  necessary.  On  the  other  hand,  the 
chordwise  distributions  of  the  bound  shroud  vortices  and  the 
kernels  will  possess  infinities.  The  coefficients  and 

will  have  square-root  singularities  at  the  leading  edge  of  the 
shroud  since  the  fluid  must  usually  flow  around  an  edge  of. 
zero  radius.  Under  special  conditions  of  tangential  flow 
entry  at  the  leading  edge,  these  singularities  will  disappear. 
The  kernels  will  always  be  singular  at  the  zero  of  their  argu- 
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ment,  however.  That  is,  in  integrating  over  the  shroud  we 
must  always  pass  over  a  contributing  element.  To  determine 
the  nature  of  these  singularities,  we  need  the  logarithmic 
behavior  of  the  Legendre  functions  near  unity.  From  Ref. 13  . 
we  have 

Q  l(^)  -  *ifn(i-l)  +  0(2-1)  fn(2-l)  (2.38) 

n-’s 

The  integrand  of  gjp»j(C,l)  ,  then,  has  a  t  £n  t  behavior 
near  the  lower  limit,  the  only  region  of  concern.  Since  this 
IS  bounded,  g  ..(die  .1)  is  regular  at  Ax..  =  C  .  But,  S...„  , 
the  remaining  terra  in  ,  gives  two  equal  logarithmic  con¬ 
tributions  and  Eq.  (?.3^)  becomes, 

VX)  -  -  i  (2.39) 

near  Ax^  C  and  the  integral  of  Eqs .  (2.35)  and  (2*36) 
must  be  interpreted  in  the  Cauchy  principal  value  sense.  This 
singularity  is  due  to  the  shed  vortices  which,  in  the  neigh¬ 
borhood  of  Ax^  -  C  ,  appear  as  a  semi-infinite  sheet  of 
vortices  of  constant  strength.  In  contrast,  differentiating 
Eq.  (2-39).  we  have  the  usual  Cauchy  singularity  of  two- 
dimensional  airfoil  theory.  Because  )/47t  is  the 

bound  shroud  vorticity  contribution,  this  result  le  antici¬ 
pated.  In  fact  with  neither  swirl  nor  propeller  terms, 

Eqs.  (2.j9)  and  (2.3l)  reduce  to  the  zeroth  harmonic  for  B 
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which  represents  an  axisynnetrlc  ring  nt  zero  angle  of  attack 

In  a  iinifora  stream  U  .  For  the  limiting  case  of  infinite 

duct  radius,  the  regular  part  of  the  kernel  vemishes,  leaving 

the  singular  p6urt  as  the  proper  result. 

Fotirth,  after  removing  the  swirl,  the  decoupled  equations 

are  identified  with  the  general  equations  for  eui  infinitely 

9 

thin  ring  wing  as  developed  by  Welssinger  .  For  the  pro¬ 
peller  and  camber  terms,  this  simply  requires  their  reinter¬ 
pretation  as  the  coefficients  of  a  trigonometric  expzuision  for 
any  specified  radial  wash  that  the  shroud  system  must  cancel. 

For  the  kernel,  though,  we  need  to  show  that  the  sum  of  the 
contributions  from  Bqs.  (2.9)  and  (2.37)  equals  his  results. 

The  Cauchy  aingul2u:ity,  as  we  have  seen,  coincides  with  the 
singularity  of  his  kernel.  This  comes  from  the  “bound"  term 
-Ax  s\/^^  •  which  is  invariant  with  the  swirl  and  which  is 

V  ItlN 

I"? 

equated  to  its  counterpart  simply  by  the  relation  between 

17 

Q^_jj  emd  the  Riegels  function  , 

(2.4o) 

where  U  =  2/(<Uj^+l)  .  To  equate  -g^(AXy,  1  )/4Tr  and  its 
counterpaurt ,  or  the  "hihed"  terms,  is  more  difficult.  One  way 
is  to  re-express  T^^  in  the  Legendre  integral  form  of 
Eq.  (2.8)  and  interchange  the  order  of  integration.  Integra¬ 
tion  on  T  represents  a  semi-infinite  straight  vortex  and 
cam  be  caurried  out.  The  remaining  Integration  then  agrees 
with  the  corresponding  v/ci  ssinger  terms  in  integral  form.  Since,  in 


turn,  he  expresses  these  in  Riegels  fvmctions,  the  final 
identity  by  means  of  Bq.  (2.40)  becomes  for  mN  0  , 


4c 


T*  dT  = 

2^ 


7T  + 


Ax  T  „ 
V  mN 


(2.41) 


a  result  which  is  needed  later.  The  arguments  on  the  RHS 
are  cSj^  . 

2.6  Shroud  and  Propeller  Loading 

After  the  determination  of  the  shroud  vortex  distribu¬ 
tions,  the  total  perturbation  velocities  induced  by  the  ducted 
propeller  system  may  be  computed  at  an  arbitrary  field  point. 
This  permits  the  calculation  of  the  shroud  pressure  loadings 
as  well  as  the  propeller  inflow  velocities. 

The  required  velocities  gp  ,  qp.  ,  and  .  are 

given  in  general  by  Eqs.  (1.12),  (1.13),  (1-1^)  and  (1.15). 

Any  component  of  the  flow  is  evaluated  by  taking  the  appro¬ 
priate  scalar  dot  product  into  a  unit  vector  in  the  desired 
component  direction.  Reduction  tc  harmonic  expressions  in 
terms  cf  the  tabulated  Legendre  functions  of  second  kind  and 
half  order  follows  the  previous  derivations. 

To  find  the  pressure,  we  use  a  linearized  form  of 
Bernoulli's  equation.  In  the  propeller-fixed  coordinates, 
the  flow  18  steady  but  rotational.  Therefore  the  total  head. 
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or  static  plus  dynamic  pressure,  is  constant  only  along  a 
streamline  .  That  is, 

p  +  pq  •  3/2\^  = 

P«  +  P(U^  +  n^r^)/2  (2.42) 

where  p  is  the  static  pressure  and  p  is  the  fluid  density. 

By  ovu'  assumption  of  small  disturbances,  the  deflection  of  any 
streamline  is  negligible  in  Eq.  (2.42),  and  we  may  use  the 
same  radius  r  for  p  and  £  as  at  infinity.  Omitting  the 
higher  order  velocity  terms  as  before,  we  have  for  dp  =  (p-p^^)  , 

dp  =  -  P£y  •  {  qp  +  3r'  Sly  Sy'  ^  ^2*^3) 

with  from  Eq.  (1.1)  for  a  =  0  .  Only  the  axial  and  tan¬ 
gential  components  of  ...  are  necessary  because  £y  does 

not  have  a  radial  component. 

From  Eq.  (2.43)  and  the  induced  velocities,  the  inner  zmd 
outer  shroud  pressure  distributions  may  be  computed.  In  derail, 
these  calculations  are  quite  complex,  but  the  net  loading,  or 
the  inner  pressure  minus  the  outer  pressure  dp]  is  simply. 


dp]  =  -pU7 


(2.44) 
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This  we  find  from  the  Kutta-Joukowski  force  on  a  bouiid  vortex 
element  or  from  Eq.  (2.43)  since  the  only  discontinuous  velo¬ 
city  component  at  the  element  is  produced  by  the  element  itself 
and  is  in  a  plane  normal  to  the  element.  The  leading  edge 
force,  if  present,  may  be  calculated  in  the  usual  manner 
from  the  coefficient  of  the  singular  term  of  the  Glauert 

series  for  A._  and  B 
m  m 

Calculation  of  the  axial  and  tangential  components  of 
the  induced  velocities  at  the  propeller  plane  will  provide 
the  inflow  quantities  corresponding  to  the  given  loading. 

With  these  quantities,  the  blade  angle  settings  and  forces 

l8 

can  be  computed  from  the  usual  propeller  blade  element  theory 
of  S,  Drzewiecki.  Since  the  inflow  is  due  to  all  the  vortex 
elements,  the  result  is  analogous  to  the  Prandtl  lifting  line 
approximation  for  wings  of  moderate  aspect  ratio  and  is  con¬ 
sistent  with  our  representation  of  the  blade  as  a  vortex  spike. 
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SOLOTIONS  FOR  DUCT  LOADING 


3.1  Fundaiaental  Solution 

The  zeroth  harmonic  of  the  shroud  vorticity  distribution 
in  propeller-fi'^ed  coordinates  is  independent  of  angular  posi¬ 
tion  on  the  shroud  for  a  given  chordwise  position.  On  the 
other  hand,  if  we  transform  to  duct-fixed  coordinates,  we  see 
that  this  harmonic  corresponds  to  the  time-independent  part  of 
the  shroud  vortex  distribution.  As  such,  it  gives  the  average 
or  fundamental  duct  loading. 

Setting  m  equal  to  zero  in  Eg.  (2.35).  then,  we  obtain. 


2€(Xg) 


Bq  {\) 


(3.1) 


where  the  kernel  ~k1  , 

o 


V  ^ 


(1  +  £^^/2) 


is  simplified  by  the  reduction  of  the  "bound"  terms  from 
Eq.  (2.22)  and  by  the  absence  of  the  "shed"  terms.  In  fact, 
Eqs.  (1.22)  cind  (2.8)  show  that  the  radial  velocity 

Influence  function  of  a  ring  vortex  of  constant  unit  strength 
evaluated  at  the  same  1  lius  R  . 

From  Eq.  (2.39),  the  singular  part  of  the  kernel  is 
l/2rrAx^  .  Removing  this,  we  hive  plotted  the  regular  part  of 
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~27TKo  ,  which  is  antisynsnetric,  in  Fig.  3*1  and  tabulated  the 

values  in  Table  3 *3-  for  C  ^  Ax^  ^  1  .  The  regular  part 

represents  the  three-dimensional  effect  of  the  vortex  ring. 

From  Eq.  (2.38) »  it  behaves  as  Ax^  £n  Ax^  at  the  origin, 

increases  to  a  maximum,  amd  then  decays  monotonically .  The 

computational  procedures  followed  those  developed  for  the 
15 

Legendre  functions  ,  and  the  results  agree  with  Weissinger's 

Q 

alternate  calculation  ,  or  •  where  t]  corresponds  to 

^  • 

The  propeller  contribut.lon  is  also  simplified  for  =  0  . 
Frcm  Eqs.  (2.22),  (2.29)  and  (2,32).  we  find 

A  priori,  we  might  not  anticipate  thic  non-periodic  contribu¬ 
tion,  but  Its  identity  will  be  established  later. 

Since  Bp  'q  appears  in  parallel  with  the  camber  terra,  we 
see  that  it  is  of  cons  ierable  significance.  Without  the  pro¬ 
peller  terra,  Eq.  (3*i)  is  the  equation  for  a  ring  wing  at  zero 
angle  of  attack,  consequently,  the  fundamental  solution  for 
the  ducted  propeller  may  be  taken  from  this  ring  wing  solution 
by  using  an  effective  camber  , 

^e  ~  ^ 


(3.3) 


^5 


-  Present  Work 

12 

- - Welssinger 


FIGURE  3.1 

BEHAVIOR  OF  REGULAR  PART  OP  KERNEL  FOR  2EROTH  HARMONIC 

iIq  =  t  1/2 


TABLE  3.1 

REGULAK  PART  OF  KERNEL  FOR  2ER0TH  HARMONIC 


-2irK'  -  \/i 

0 

o.co 

o.ocooo 

C.05 

0.07981 

C.IO 

0.1329c 

C.15 

0.17629 

0.20 

0.21305 

0.25 

0.24482 

0.30 

0. 27253 

C.35 

0.29683 

0.40 

0.31815 

0.45 

0.33684 

0.50 

0.35322 

0.55 

0.36751 

0.60 

c. 37989 

0.65 

C. 39058 

0.70 

0.39966 

0.75 

0.40734 

0.80 

0.41370 

0.85 

0.41889 

0.90 

0.42298 

0.95 

0.42610 

1.00 

0.42833 
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In  other  words: 

The  me^m  pressure  difference  across  the  shroud  of  a  ducted 
propeller  at  zero  incidence  in  a  uniform  flow  is  identical  to 
that  of  a  similar  ring  wing  of  different  canfcer. 

Thus,  experimental  measurements  on  ring  wings  by 

iq 

0.  Ladurner  and  others,  the  rheoelectrlc  analogy  of 
L.  Malavard^*^,  and  the  theoretical  methods  of  Weiss Inger^ 
can  be  applied  to  ducted  propellers.  Two  points  should  be 
noted.  First,  this  result  is  exactly  valid  only  within  the 
approximations  of  a  linearized  theory  and  of  an  invlscid 
fluid.  Second,  the  shroud  does  not  have  to  be  of  zero 
thickness.  Though  they  do  alter  the  effective  candser,  thick¬ 
ness  effects  can  be  superimposed^^ within  the  limitations 
of  thin  airfoil  theory.  In  fact,  with  respect  to  the  shroud 
thickness,  they  are  Identical  for  a  ducted  propeller  and  its 
equivalent  ring  wing. 

3-2  Effect  of  Parameters  on  the  Fundamental  Load 

The  modification  to  the  effective  camber  by  the  propeller 
for  the  fundamental  solution  depends  on  the  following  parameters 

( 1 )  Blade  number 
(li)  Blade  circulation  strength 

(iii)  Advance  raMo 

(iv)  Tip  clearance 

(v)  Propeller  position 

(vl)  Strength  of  the  vortices  shed  from  the  blades 
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To  the  extent,  then,  that  camber  effects  are  known  for  the 
ring  wing,  the  effects  of  these  parameters  on  determine 

their  influence  on  the  steady  shroud  load. 

To  study  we  express  the  blade  circulation  in  a 

general  form,  in  particular,  in  most  cases  for  no  loading  at 
the  root  and  finite  tip  clearemce,  the  circulation  distribu¬ 
tion  will  be  zero  at  r^  =  C  ,  without  a  hub,  and  at  r^  =  H  . 
Therefore,  we  take 


r  = 


sin  JT  r 


(3.4) 


where  ,  for  convenience,  is  the  maximum  value  of  F  and 
r^  =  r^/M.  .  From  Fourier  emalysis,  we  have, 


F 

J 


rT-  sin  J7T  r  .  dr 

■’  V  V 


(3.5) 


to  determine  the  coefficients. 

With  the  blade  circulation  distribution  of  Eg.  (3.4)  the 
propeller  contribution  to  the  effective  camber  becomes 


®F'0 


J 


Xj(^"p.4) 


(3.6) 


after  interchange  of  the  order  of  summation  cind  integration. 
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The  chauracterlstlc  ftinctions  Xj  given  by 


=  1  + 


_  o  ,  ^  v2 

Ax  +  (l-iir 

■  kl.m  . 

2ur„ 


(3.7) 


In  the  form  of  Eq.  (3*6),  the  dependence  of  on  the 

various  parameters  can  be  examined. 

The  simplest  dependence  is  on  the  blade  nupaber,  magni¬ 
tude  of  the  blade  loading  and  the  advance  ratio.  For  the 
first  t%ra,  the  propeller  camber  contribution  varies  linearly 
in  a  uniform  fashion  along  the  duct.  Nhen  N  or  is 

zero,  Bp.Q  is  zero  as  required.  On  the  other  hand,  it  is 
Inversely  proportional  to  the  advance  ratio.  As  J 0 
either  by  U -♦  C  or  n-**»  ,  Bp'Q  becomes  Infinite,  which 
corresponds  to  the  collapse  of  the  shed  propeller  vortices 
into  a  disk  of  concentric  vortices  of  infinite  strength  at 
the  propeller  plane.  In  contrast,  when  the  advance  ratio 
becomes  infinite,  the  caudaer  modification  disappears,  we  may 
interpret  this  result  in  the  following  way.  Bp'o  depends  only 
on  the  variation  of  the  blade  circulation  and  not  on  the  blade 
circulation  itself.  That  is,  the  only  reason  we  have  an  effect 
at  all  for  finite  advance  ratios  is  the  essential  "cJtearing" 
of  the  basic  periodic  character  of  the  blades  by  their  shed 
vortices  as  they  rotate  out  of  the  radial -axial  planes  in 
which  they  were  shed.  For  infinite  advamce  ratio,  the  shed 
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vortices  trail  back  aociaily.  Then  there  is  no  smeauring  amd 
consequently  no  effective  camber. 

The  Influence  of  the  three  r^aining  parameters,  the 
blade  circulation  profile  or  r^'s  ,  the  propeller  position 
3^  ,  and  the  tip  clearance  variable  p.  is  not  so  easily 
obtained.  Analytically,  we  see  from  Eq.  (3.7)  that  each 
characteristic  function  is  sysmetric  about  the  propeller 
plane.  Since  the  corresponding  Integration  of  Xj  from 
0  to  icg  is  cintlsyiiinetric,  the  equivalent  shroud  shape 
"induced"  by  the  Fourier  component  of  the  blade  load¬ 
ing  is  antisymmetric  with  an  inflection  point  at  =  0  . 

Differentiation  of  Xj  with  respect  to  Ax^  shows  that  the 
slope  of  this  induced  casiberllne  or  contribution  to  the  effec¬ 
tive  camber  is  a  finite  maximum  or  minimum  at  the  propeller 
plane;  i.e.,  dS^/dAXp  =  AXp/ur^  -*  0  .  The  result  is  valid 
for  the  complete  ramge  of  0  ^  U  ^  1  as  verified  by  the 
logarithmic  behavior  of  the  Legendre  function  near  its  singu¬ 
larity  and  its  asymptotic  behavior^^  of  ••  0  • 

This  agrees  wirh  the  usual  streamline  patterns  depicted  for 
a  propeller  which  show  larger  radial  velocities  nearer  the 
propeller. 

To  furnish  further  insight  as  well  as  practical  %irorking 
figures,  we  have  computed  the  characteristic  functions  for 
j  =  1,2, 3,^,5  at  values  of  u-  =  C.75  .  0.90  ,  C.95  and  C.99 
If  the  propeller  plane  is  kept  within  the  shroud,  the  maximum 
range  of  AXp  required  for  any  calculation  is  2X  .  We  have 
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considered  values  of  C  ^  or  a  minimum  of  X  =  Oo 

for  the  case  of  the  propeller  plane  at  the  shroud  entreUice 
or  exit.  For  a  centrally  located  propeller,  the  m2ucimura  of 
X  =  1  is  attained. 

The  results  are  given  in  Figs.  3.2,  3-3.  3.^.  and  3.5 
along  with  their  tabulated  values  in  Tables  3.2,  3.3,  3.4  and 
3.5.  AS  previously  noted,  the  greatest  effects  occur  at  the 
propeller  plzme  and  decay  with  Increasing  axial  separation. 

The  decay  is  faster  for  the  higher  harmonics.  Analogous 
behavior  is  found  for  m-  or  radial  separation  and  is  illus¬ 
trated  in  Fig.  3-6,  a  cross  plot  of  the  absolute  values  of 

_  _ii 

Xj  at  AXp  =  C  .  The  tabulation  is  accurate  to  10  or  better. 

As  the  tip  clearance  decreases  we  know,  for  a  constant 
geometry,  that  the  load  near  the  tip  does  not  generally  decay 

23 

so  early  and  that  the  subsequent  shed  vortex  grows  in  strength  . 
In  the  ideal  flow,  however,  this  vortex  disappears  when  the 
tip  just  touches  the  shroud  and  r'=  0  .  The  effect  of  the 
tip  vortex,  then,  may  be  examined  by  assuming  a  constant  blade 
circulation  throughout  the  blade  and  approximating  the  tip 
beh'vior  by  a  step  function.  Substituting  f  = 
into  Eq.  (3-2)  and  evaluating  th2  integral,  we  get 


TO 


2 

2-n  J 


1  K 


2u 


(3.8) 


consequently,  the  contribution  of  a  tip  vortex  of  strength 

r,,  to  the  effective  shroud  camber  ci?-.  be  found  from  the 
H 
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TABLE  3.2 

CBAEACTBRX8TIC  FQKTICBn  Xj  TOE  |i  «  0.75 


Xg 

0.00 

0.25560 

0.08093 

0.05 

0.25317 

0.07963 

0.10 

0.24639 

0.07602 

0.15 

0.23614 

0.07077 

0.20 

0.22356 

0.06464 

0.25 

0.20972 

0.05831 

0.30 

0.19543 

0.05219 

0.35 

0.18126 

0.04654 

0.40 

0.16758 

0.04147 

0.45 

0.15462 

0.03700 

0.50 

0.14248 

0.03308 

0.55 

0.13121 

0.02968 

0.60 

0.12081 

0.02673 

0.65 

0.11126 

0.02416 

0.70 

0.10249 

0.02192 

0.75 

0.09447 

0.01996 

0.8o 

0.08714 

c. 01824 

0.85 

0.06044 

0.01673 

0.90 

0.07433 

0.01538 

0.95 

0.06874 

0.01418 

1.00 

0.06363 

0.01511 

0.03917 

0.02296 

0.01503 

0.03837 

0.02242 

0.01466 

0.03619 

0.02098 

0.01364 

0.03309 

0.01897 

0.01224 

0,02962 

0.01678 

0.01075 

0.02617 

0.01467 

0.00934 

0.02299 

0.01277 

0.00810 

0.02018 

0.01114 

0.00705 

0.01775 

0.00976 

0.00617 

0.01569 

0.00861 

0.00544 

0.01395 

0.00765 

0.00483 

0.01247 

0.00684 

0.00433 

0.01121 

0.00616 

0.00390 

0.01013 

0.00557 

0.00353 

0.00920 

0.00507 

0.00321 

0.00840 

0.00463 

0.00294 

0.00769 

0.00425 

0.00270 

0.00707 

0.00391 

0.00248 

0.00652 

0.00361 

0.00229 

0.00602 

0.00334 

0.00212 

0.00558 

0.00310 

0.00197 
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TABLE  3.3 

CBAXACTESISTIC  PUMCTIONS  Xj  J’OR  H  =  0.90 


'JL 

£2. 

0.00 

0.39179 

0.14562 

0.05 

0.38341 

0.13947 

0.1c 

0.36257 

0.12492 

0.15 

0.33623 

0.10797 

0.20 

0.30873 

0.09194 

0.25 

0.28215 

0.07800 

C.30 

0.25734 

0.06632 

0.35 

0.23461 

0.05671 

0.40 

0.21396 

0.04885 

0.45 

0.19529 

0.04243 

0.50 

0.17845 

0.03716 

0.55 

0.16328 

0.03282 

0.60 

0.14960 

0.02921 

0.65 

0.13727 

0.02620 

C.70 

0.12612 

0.02365 

0.75 

0.11605 

0.02149 

0.80 

0.10693 

0.01962 

0.85 

0.09865 

0.01800 

o.yo 

0.09114 

c. 01659 

0.95 

0.08431 

0.01535 

l.CC 

0.07809 

0.01424 

0.07773 

0.04866 

0.03338 

G.073C5 

0.04499 

0.03044 

0.06246 

0.03701 

0.02426 

0.05100 

0.02892 

0.01834 

C. 04107 

0.02241 

0.01386 

0.03318 

0.01760 

0.01073 

0.02713 

0.01415 

0.00859 

0.02254 

0.01167 

0.00710 

0.01906 

0.00988 

0.00604 

0.01639 

0.00853 

0.00525 

0.01431 

0.00751 

0.00465 

0.01266 

0.00670 

0.00417 

0.01133 

0.00604 

0.00378 

0.01024 

0.00550 

0.00345 

0.00933 

0.00504 

0.00317 

O.OO855 

0.00464 

0.00293 

0.00788 

0.00430 

0.00272 

0.00729 

0.00399 

0.00253 

0.00678 

0.CC372 

0.00236 

0.00631 

0.00347 

0.00220 

c. 00589 

0.00325 

0.00206 
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TABLE  3.^ 

CHARACTERISTIC  FUNCTIONS  Xj  FOR  U  =  0-95 


+Ax„ 

r-,.B 

III 

is. 

0,00 

0.46299 

0.18750 

0.05 

0.44552 

0.17310 

0.10 

0.4105C 

0.5.4641 

0.15 

0.37306 

0.12091 

0.20 

0.33766 

0.09956 

C.25 

0.30543 

0.08238 

C.30 

0.27652 

O.C6874 

0.35 

0.25071 

c. 05795 

0.^0 

0.22773 

0.04939 

0.45 

c. 20725 

0.04256 

0.50 

0.18898 

0.037C7 

0.55 

0.17264 

0.03262 

C.6C 

O.158CI 

0.02898 

0.65 

0.14488 

0.02596 

0.70 

0.13307 

0.02344 

0.75 

0.12241 

c .02131 

0.80 

0.11279 

0.01948 

0.85 

0.10408 

0.01791 

0.90 

0.09618 

0.01654 

0.95 

0.08900 

0.01533 

1.00 

0.08247 

0.01426 

il 

zh. 

0.1C654 

0.07000 

0.04994 

0.09451 

0.05980 

0.04120 

0.07 381 

0.04347 

0.0281c 

0.05605 

0.03078 

0.01882 

0.04279 

0.02225 

0.01314 

0.03326 

0.01669 

0.00972 

0.02647 

0.01306 

0.00762 

0.02161 

0.01064 

0.00628 

0.01809 

0.00898 

0.00537 

0.015^^8 

0.00779 

0.00472 

0.01351 

0.00690 

0.00423 

0.01198 

c. 00621 

0.00384 

0.01076 

0.00565 

0.00352 

0.00377 

0.00519 

0.00325 

0.00895 

0.00480 

0.00302 

0.00825 

0.00445 

0.00281 

O.OOf65 

0.00415 

0.00262 

0.00712 

0.00388 

0.00245 

0.00-664 

0.00364 

0.0023c 

0. 00622 

0.00341 

0.00216 

0.00583 

0.00321 

C. 00204 
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TABLE  3.5 

CHARACTERISTIC  FUNCTIONS  X  FOR  M-  ■=  0.99 


2k 

0.00 

0.55041 

0.24839 

C.05 

0.50412 

0.20640 

C.IC 

0.44968 

C •16?77 

C.15 

C. 40135 

0. 12906 

C.2C 

c. 3592c 

0.10346 

0.25 

0.32248 

0.08400 

C.30 

C. 2904 3 

0. 06913 

C.35 

c. 26236 

0. 05769 

C.4G 

o.2j^68 

0.04981 

C.45 

0.21590 

0.04185 

C.5C 

0.19661 

0.03633 

C.55 

0. 17945 

0. 03191 

c  .60 

0.16415 

0.02832 

C.65 

C.1504D 

0.02538 

C.fC 

0.13818 

0.02293 

C.75 

0.12712 

0 . 02088 

C.8C 

G. 11715 

0.01913 

C.85 

C.  10813 

0.01762 

c  .9c 

0.09996 

G.VI631 

0.09253 

0.01516 

l.CC 

C. 08578 

0.01413 

0.15401 

0.10893 

0.08285 

0.11589 

0.07423 

0.05118 

0.08094 

0.04614 

0.02850 

0.05748 

0.02973 

0.01693 

0.04198 

0.02025 

0.01105 

0. 03166 

0.01468 

0.00796 

0.02471 

0.01132 

0.00625 

C.CI995 

0.00921 

C.CO523 

0.01662 

0.00783 

0.00457 

0.01423 

C.OO6S9 

0.00412 

0.01246 

0.00619 

O.CO376 

O.Cllll 

0.00565 

0.00348 

O.OICO5 

C.CC521 

0.00324 

0.00920 

0.00484 

0. 00303 

0.00848 

0.00452 

0.00284 

0.00788 

0.00423 

C.CC267 

0.00735 

0.00398 

0. 00251 

0.00688 

0.00374 

0.00237 

0.00646 

0.00:553 

0.00224 

0.00607 

0.00333 

G. 00211 

C.(.0572 

0.00315 

O.C020C 
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tabulated  values  of  tiie  Legendre  function^^' .  Over  a 
range  of  values  of  u  and  ,  the  results  may  be  esti¬ 

mated  from  Fig.  2.1.  For  zero  tip  clearance,  they  possess 
a  logarithmic  singularity  near  the  propeller  plane.  Since 
Eqs.  (1.6),  (1.18)  and  (2.S)  show  that  the  radial  velocity 
at  ,  R  induced  by  a  ring  vortex  of  radius  r  is 

e  Y 

.  a  perfect  differential  in  , 

®q.  {3»8)  ia  equivalent  to  the  effect  of  a  semi-infinite 
cylindrical  distribution  of  vortex  rings  of  radius  S  , 

Jr 

extending  from  x^  to  +“  .  The  corresponding  constant 
strength  per  unit  length  in  the  axial  direction  is  UNrj^u/xJ  , 
or  alternatively,  a  uniform  slipstream  velocity  increment  of 
the  same  amount.  In  this  way,  the  previous  results^  using 
the  method  of  singularities,  i.e.,  determing  the  effect  of 
the  propeller  by  a  pressure  Jump  and  uniform  slipstream,  can 
be  related  to  the  present  analysis  by  setting  u  =  1  and 
Xp  *  X  in  Eq.  (3 *8). 

3.3  Limiting  Case  of  Infinite  Advance  Ratio 

As  the  forward  velocity  U  becomes  very  large  compared 
to  the  rotational  velocity  HRp  at  the  propeller  tip,  we  have 
shown  that  the  swirl  terms  disappear  from  the  kernel.  This 
decoupled  cur  governing  equations,  Eqs.  (2.35)  and  (2.36), 
for  the  chordwlse  shroud  vortex  distributions  and  reduced  the 
kernel  Kjjj  to  Weissinger’s  result.  In  terms  of  Legendre  func¬ 
tions,  the  reduced  kernel  is  given  by  the  sum  of  Eqs.  (2.9) 
and  (2.37). 
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NOW,  if  we  let  the  advance  ratio  become  infinite  in  the 
"shed"  propeller  terra  Cp^  and  follow  the  derivation  of 
Eq.  (2.37).  we  find  that  this  term  is  purely  imaginary,  or 


The  argument  of  T^^  is  equal  to  with  Jp.?  replaced  by 

its  equivalent  t  .  From  their  definitions,  then,  is 

identically  zero  in  this  case  and  hp^  is  given  as  minus 
twice  the  imaginary  part  of  Eq.  (5.9).  The  single  "bound" 
propeller  terra  remains  unchemged. 

With  these  results,  Eqs.  (2.35)  Jind  (2.36)  are  fxirther 
simplified  for  an  infinite  advance  ratio.  The  subsequent 
shroud  vortex  distributions  for  any  harmonic  can  be  found 
from  Weissinger’s  solutions^  for  a  ring  wing  at  zero  angle 
of  attacl;  in  a  uniform  stream.  In  particular,  for  we 

see  that  —  are  all  zero  since  the  real  shroud 

f'amb'r  appears  only  in  the  zeroth  harmonic.  The  remaining 
B^  with  Bj-,r  -  C  is  the  distribution  for  the  act.al  duct, 
which  is  assumed  axisyimnetric ,  without  the  propeller  present. 
This,  of  course,  is  the  same  as  noted  for  B^  for  finite  J 
since  there  are  no  shed  vortiocs  for  axisyraraetric  loading. 

On  the  other  hand,  for  A  ,  A^  =  0  and  A, ,A^,  ... 

rn  o  X  <2 
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I 

i 

I 


correspond  to  the  vortex  distributions  of  cin  equivalent 
asymmetric  duct  under  identical  flow  conditions.  The  equiva¬ 
lent  axially  asymmetric  camber,  say  e  ,  has  just  a  sine 

u  ** 

harmonic  dependence,  or, 

00 

^jm  (^p)  (3.10) 

m=l 

where  the  axial  candaer  coefficients 


€ 


jm 


=  fcim 
j  -*  ® 


(3-11) 


follow  directly  from  Eqs.  (2.30)  and  (3.9).  That  is: 

At  high  advance  ratio,  the  resultant  shroud  loading  on 
a  ducted  propeller  at  zero  emqle  of  attack  in  a  uniform 
inviscid  flow  is  equal  to  '-ha  load  on  the  isolated  shroud 
Plus  the  loading  on  an  equivalent  asymmetric  ring  wing. 

The  inner  integration  on  T^(a)^)  over  t  in  Eq.  (3.9) 
is  given  by  Eq.  (2.41)  for  r^  =  1  Juid  by  the  infinite  series. 


"  Aii 


(J.l)  . 


•nr 


£  I  (f^-f-"^)(£-^-f) 


(®c 


2d 


2, mu 


(3.12) 
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for  0  <  1  .  The  argvuaent  of  the  (j+1)^^  derivative 

is  5^  =  {l+(l-g^)‘=^/2j^]  or  of  Eq.  (2.18)  with  Ax^  =  0  , 

i.e.,  the  propeller  plane,  zmd  f  =  {a)^-(S^  -1)*]  .  To  derive 
(3*12) ,  Tjjijj  is  re-expressed  in  the  integral  form  of  the 
Legendre  function  of  Eq.  (2.8)  emd  the  order  of  integration 
is  interchamged .  The  integration  on  t  ,  \diich  represents  the 
influence  of  a  straight  semi-infinite  vortex,  can  then  be 
performed.  Expanding  part  of  the  subsequent  integrand  over 
9  in  reciprocal  square-root  powers  of  {(l-r^)^+4r^,  sin^'’9)  , 
the  integral  of  each  term  is  identified  from  Eq.  (2.8)  as  a 
derivative  of  a  Legendre  function  of  second  kind  and  half 
order.  The  other  part  is  evaluated  by  suitable  contour  inte¬ 
gration. 

3.4  Infinite  Blade  Nunber 

In  rnamy  instances,  the  propeller  blade  nunber  may  become 
quite  large.  For  this  case,  an  actuator  disk  solution  is 
desired  which  is  the  proper  limit  of  our  equations,  retaining 
the  radial  propeller  load  variation. 

The  propeller  disk  loading  is  proportional  to  the  product 
of  the  blade  nunber  and  the  strength  of  tl»e  blade  circulation. 
TO  obtain  an  actuator  disk  solution,  we  let  the  blade  nunber 
approach  infinity  but  decrease  the  strength  of  the  circulation 
such  that  the  product  remains  finite,  or 


M 


t.im  ur 
M  -  a  ,  0 


(3-13) 


If  we  examine  the  propeller  modification  of  the  shroud  camber 
by  the  zeroth  harmonic,  see  Eq.  (3*6),  we  find  that  for  these 
parameters  its  only  dependence  is  linear  with  their  product. 

In  the  limit,  then,  this  camber  modification  persists  unchanged 
and  equal  to  the  finite  blade  result. 

On  the  other  hand,  the  frequency  of  the  first  harmonic 
and  other  time  dependent  terms  of  the  shroud  loading  becomes 
infinite.  By  definition  of  any  actuator  disk  model,  however, 
the  flow  field  must  be  time  independent  in  coordinates  trans¬ 
lating  with  the  propeller.  Since  these  higher  harmonics  are 
purely  sinusoidal,  their  amplitude  fcr  infinite  blade  number 
and  finite  disk  loading  must  approach  zero.  Physically,  this 
establishes  that: 

The  average  pressure  distribution  on  the  shroud  of  a 
ducted  propeller  with  a  finite  blade  number  corresponds  to 
the  total  loading  on  a  similar  model  with  an  infinite  blade 
number  but  the  same  disk  loading. 

For  an  analytic  proof,  we  have  to  note  that  the  parameters 
ra,N,  and  occur  in  the  propeller  terms,  see  Eqs.  (2.29) 
(2.3c)  and  (2.31),  only  as  the  product  mK  and  .  The 

kernels,  of  course,  contain  just  mW  as  a  parameter.  liathe- 
raatically,  then,  we  have  to  expand  for  large  order  near 

its  singularity  to  show  that  all  the  propeller  terras  except 
Bpj/  vanish  and  that  the  integral  equations  decouple. 

The  zeroth  harmonic  of  the  propeller  terms  for  the  actu^. - 


tor  disk  limit  has  a  broader  significance,  ihe  term 


66 


gives  the  radial  velocity  component  of  a  generalized  actualcr 
disk.  The  computation  of  the  axial  and  tcmgcntial  velocity 
components  is  analogous  to  that  of  .  This  model  is  not 

limited  to  a  uniform  disk  loading,  nor  to  zero  angle  of  attack 
if  the  slipstreeim  deflection  is  known.  In  other  v/ords; 

The  steady  part  of  the  flow  field  of  any  propeller  with 
finite  blade  number  corresponds  to  the  generalized  actuator 
disk  solution  with  the  same  disk  loading . 

3*5  General  solution  for  Higher  Harmonics 

we  have  shown  that  the  higher  harmonics,  m  =.1,2,3  •••  • 
correspond  to  the  time -dependent  terms  of  the  shroud  loading. 

The  governing  equations  have,  in  general,  two  essential  fea¬ 
tures  which  make  a  nijmerical  solution  more  complex  than  the 
ring  wing  problau.  One  is  the  coupling,  and  the  other,  the 
occurrence  of  a  logarithmic  as  well  as  a  Cauchy  singularity. 

These  difficulties  can  be  treated  by  modification  of  Weissinger's 
technique,  but  further  reduction  is  more  desirable. 

First,  It  is  possible  to  remove  the  coupling  because  of 
the  "derivative"  relationship  betv/een  the  kernels.  For  m  ^  1  , 
we  let 

•“X 

=  /  B  (it  )  K  (Ax  )  dx 

B  m'  v'  m'  v'  v 


(3-L4) 
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and  note 

Sf-  V*-\>  (3-15) 

s 

With  Eqs.  (3.14)  eind  (3-15) .  Eqs.  (2.35)  and  (2.36)  become  a 
system  of  linear  differential  equations.  Or,  without  the 
shroud  camber, 


B 


rm 


01.,  ■-  I, 


(3-16) 


where  differentiation  across  the  integrals  and  ig  is 

justified  in  view  of  the  weal;  square-root  and  logarithmic 

singularities  of  the  integrands.  Multiplying  the  first  of 

Eqs.  (3  1  •}  ^  raN/Ju  and  adding  it  to  the  derivative  of 

the  second,  v;e  isolate  •.  reversing  the  procedure  and 

subtracting,  I_  .  Both  are  of  tl.e  form 
B 

l'%  o'  I  •  P(Xg)  (3.17) 


which  is  the  second-order,  linear  differential  equation  for 
a  simple  harmonic  oscillator  subject  to  a  forcing  function. 
'That  IS,  we  vie’//  the  integrals  and  oi  our  unhJiovm 

vortex  uistribuuions  as  "particle  displacement"  and  the  axial 
coordinate  as  "time".  The  "motion"  is  undait.pea  with  a 
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"natural  frequency"  a  .  The  "forcing  f mictions"  for 
and  Ig  are , 

^  “  »rr; 

’  '“ml '  -  (3-18) 

respectively. 

From  these  differential  equations.  and  Ig  are 

determined,  providing  in  turn,  decoupled  equations  for 

and  .  'fo  find  the  inhomogeneous  solution  oi  Eq.  (3.17). 

we  take  the  "indicial  response"  of  the  system  to  the  unit 

step  function  P  --  l(x  )  . 

s 

(l,c2)(l-cos  oXg)  i(rCg)  (3.19) 

24 

and  use  Duhariel's  linear  superposition  method.  After  inte¬ 
gration  by  parts,  we  have 

’  r  “3  _  _  _  _ 

J  p(:0  sin  a(Xg-x)  dx  (3.2C) 

where  P  '-.'o  '.nov^  is  v/eli -behaved  throughout  the  interval  from 
-X  to  X  for  C  <  d  <  1  . 

Eq.  (3*2^).  together  v/ith  Eqs .  (3.14)  and  (3.I0),  yield 
a  simpler  system  of  decoupled  equations  for  the  time-dependent 
shroud  vortex  distributions  v.’hich  correspond  to  the  first  and 


higher  harmonics  of  the  radial  flow  induced  by  the  propeller. 
The  form  of  these  equations  is  similar  to  the  general  ring 
wing  problem,  but  the  kernel  _s  completely  different.  However, 
if  we  differentiate  with  respect  to  ,  carry  the  diffcren- 
tiation  back  inside  euid  ,  and  use  Eq.  {Z-lj),  we  find 

for  the  inhomogeneous  solution. 


(x)  COS  0  Ax  dx 
A 


(3.21) 


-X 


and, 


)  cos  a  hit  dx 
E'  ' 


I 

-X 


B  (x  )  K  (Ai:  )  die 
m'  v'  m'  V  V 


(3.22) 


v/herc  the  "bound"  term  c£  k'  ,  or  -S ,  is  the  same  as 

m  mH 

for  the  ring  wing.  Numerical  solutions,  then,  can  follow 

9 

c.ssentially  along  the  lines  of  weissinger's  . 

The  decoupling  technique  which  we  have  employed  raises 


an  import£mt  question.  In  addition  to  the  inhomogeneous 
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solution  of  the  differential  equation  of  Eq.  (3.17),  there 
is  the  homogeneous  solution  , 

cos  a  Ax  I  + 

(l/o)  sinaAxj  'x3»23) 

where  denotes  the  value  of  any  reference 

point  between  the  leading  and  trailing  edges  with  the  leading 
edge  excluded;  and  Ax!  ,  the  distance  from  th.-  reference 
point  to  Xg  .  These  solutions  introduce  a  total  of  four 
unknown  consteints  for  the  corresponding  I.  and  I  .  Two 
of  these,  though,  result  from  the  differentiation  required 
for  decoupling.  They  may  be  removed  by  resubstitution  into 
eitl'.or  of  Eqs.  (3*16)  .  -...ch  shov/s  that  the  relations, 

must  be  fulfilled  in  order  for  the  coupled,  homogeneous  equa¬ 
tions  to  be  satisfied.  We  are  thus  left  with  tv/o  undetermined 
constants  as  we  should  expect. 

That  is,  consider  an  a::isvnimetric  ring  v;inq  at  zero  inci¬ 
dence  in  a  uniform,  translating  and  rotating  flov/,  with  the 
axis  or  rotation  coincident  With  the  wing  axis .  The  equations 
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8««n  to  Indicate,  then,  that  It  is  possible  in  an  inviscid 
fluid  to  have  self-sustaining  load  distributions  of  periodic 
nature.  Or,  the  radial  velocity  induced  by  the  helical 
vortices  shed  frosi  the  shroud  for  certain  shroud  vortex  dis¬ 
tributions  balances  the  velocity  induced  by  the  bound  shroud 
vortices,  litis  is  independent  of  tJie  casiber  within  the 
present  assunptions.  but  dependent  on  the  swirl,  vanishing 
as  the  shroud  "advance  ratio"  becosies  infinite  and  decouples 
tlte  governing  equations.  The  possible  existence  of  these 
self-sustaining  solutions  appears  to  be  similar  in  nature  to 
the  .’ece.it  investigation  of  H.  Ludwieg^®,  where  he  has  ex¬ 
amined  the  ""itsbility  against  the  formation  of  helical 
vortices"  in  an  annular  flow  between  two  coaxial  cylinders. 

3«6  Example 

To  illustrate  the  practical  application  of  the  theory 
we  have  computed  the  average  shroud  losd  distribution  for 
a  particular  ducted  propeller. 

The  values  used  for  the  parameters  and  the  propeller 
loading  assuiced  are, 

J  =  1 

K  2 

Xp=0 


* 


€  =  -0.38  X 
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r  =  (2/3V3){2  Bin  Trr  -  sin  27r  r) 
k  =  0.5 


where  r  is  shovm  in  Fig.  3*7. 

With  =.  V3'/F.  ^2  =  -  2/3'/!  and  r„  =  1  .  we 

have  computed  and  plotted  in  Fig.  3»8  the  propeller  contribu¬ 
tion,  -Bp'q/2  ,  to  the  effective  camber  by  means  of  Eq.  (3*6) 
and  Table  3.4.  The  maximum  camber  distortion  which  occurs  at 
the  mid  chord  or  propeller  plane  is  about  40^5  of  the  given 
shroud  camber  at  the  trailing  edge.  The  assumed  F  corres¬ 
ponds  to  a  disk  loading,  neglecting  inflow,  of  4l.l  psf  at 
U  =  200  raph  amd  STP. 

From  Eq.  (2.44),  the  net  average  pressure  distribution 
4>p ]  becomes  by  meauis  of  Eqs.  (2.1)  and  (2.2), 


Cp]  =  -Bq(x)  (3.26) 

where  c  ]  =  Ap]/(pU"/2)  .  For  tho  determination  of  B^(x) 

^  25 

for  our  example,  we  have  employed  an  iteration  technique 

instead  of  Weissinger's  procedure  for  the  solution  of 
Eq.  (3.1)-  Tne  final  result  is  shown  in  Fig.  3 •9-  It  re¬ 
quired  eleven  non-zero  terms,  including  the  singular  term, 
of  the  Glauert  series  and  two  iterations  to  obtain  accuracy 


I 


to  three  deciiaal  places. 


BIADE  CIBCUIATION  DISTRIBUTION  FOR 


PROPELLER  CONTRIBUTION  TO  EFFECTIVE  SHROUD  CAMBER  FC» 


I 
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} 


A 


I 


7t 

Without  the  propeller,  the  net  leading  is  synnnetric;  aboit 
the  center  of  the  shroud  and  vanishes  at  the  leading  zmd 
trailing  edges  because  the  camber  is  antisymmetric.  The 
contributlcn  of  the  propeller,  however,  is  always  symmetric 
and  so  in  this  case  introduces  flow  around  the  leading  edge 
and  a  leading  edge  singularity.  This,  in  turn,  gives  rise 
to  a  leading  edge  suction  force  which  can  l^e  found  from  th' 
value  of  0.187  for  the  coefficient  of  the  singular  term  of 
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CONCLUSIONS 


A  three-dimensional  theory  has  been  developed  for  the 
diXcted  propeller  with  finite  blade  number  in  uniform  motion 
t^-.rough  em  inviscid,  incompressible  fluid  at  zero  incidence. 
Within  the  approximations  of  a  lightly-loaded  propeller  and 
of  thin  airfoil  theory,  the  following  generalizations  wmze 
obtained: 


The  effect  of  shroud  camber  appears  only  in  the 
steady  shroud  load. 

The  steady  pressure  difference  across  the  shroud 
is  identical  to  that  of  a  similar  ring  wing  of 
different  camber. 

At  high  advance  ratio,  the  resultant  shroud  loading 
is  e<|ual  to  the  load  on  the  isolated  shroud  plus 
the  loading  on  an  equivalent  asymmetric  ring  wing. 

The  steady  pressure  distribution  on  the  shroud 
corresponds  to  the  total  loading  on  a  similar 
configuration  -.'ith  infinite  blade  number  but  the 
same  radial  disk  loading. 

The  steady  part  of  the  flow  field  of  any  propeller 
with  finite  blade  number  corresponds  to  a  general¬ 
ize’  actuator  disk  solution  with  the  same  radial 
die-,  loading. 


The  theory  is  readily  applicable  to  practical  calculations, 
partic’ ; arly  the  determination  of  the  zeroth  harmonic  which 
gives  ■  he  average  shroud  loadlna.  By  variation  of  paiametern, 
the  optxmuin  configurations  of  Dlckmann  and  Weissinger^  and 
R.  B.  ^;ray  and  W.  Castles.  Jr.^^  may  be  further  investigated. 

.'lese  results  clearly  justify  similar  analytic  studies 
of  t'.e  more  complex  ducted  propeller  problems. 
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